Current noise in molecular junctions: effects of the electron-phonon interaction 
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We study inelastic effects on the electronic current noise in molecular junctions, due to the cou- 
pling between transport electrons and vibrational degrees of freedom. Using a full counting statistics 
approach based on the generalized Keldysh Green's function technique, we calculate in an unified 
manner both the mean current and the zero-frequency current noise. For multilevel junctions with 
weak electron-phonon coupling, we give analytical formulas for the lowest order inelastic correc- 
tions to the noise in terms of universal temperature- and voltage-dependent functions and junction- 
dependent prefactors, which can be evaluated microscopically, e.g. with ab-initio methodologies. 
We identify distinct terms corresponding to the mean-field contribution to noise and to the vertex 
corrections, and we show that the latter contribute substantially to the inelastic noise. Finally, we 
illustrate our results by a simple model of two electronic levels which are mutually coupled by the 
electron-phonon interaction and show that the inelastic noise spectroscopy is a sensitive diagnostic 
tool. 

PACS numbers: 72.70.-|-m, 72.10.Di, 85.65.+h, 73.63.-b 



I. INTRODUCTION 

Recent progress in experimental techniques, such as 
break junctions and scanning tunneling microscopy, al- 
lows to electrically contact single molecules, to create 
and manipulate atomic wires, and to investigate the elec- 
tronic transport properties of these nanoscopic objectsi^ 
Contrary to larger devices, atomic-scale ones usually re- 
tain their microscopic features, which are then observable 
in the transport spectroscopy. Apart from the purely 
electronic contributions, other degrees of freedom such 
as vibrational modes or local spins can be addressed and 
revealed by point-contact spectroscopj^ (PCS) or by in- 
elastic electron tunneling spectroscopy (lETS)i^ These 
techniques have been extensively used in the recent past 
to reveal inelastic features in the non-linear conductance 
due to vibrations^"— or local spin excitations,—"— trig- 
gering an intense theoretical activity. So far most stud- 
ies focused on the current-voltage characteristics, and 
present day PCS/IETS theoriesr^"— often based on ab- 
initio calculations, allow to make detailed predictions for 
the conductance that compare favorably with experimen- 
tal results. 



Electronic current (shot) nois( 



■ 32.33 



is another quan- 



tity of fundamental interest in transport, representing 
the second cumulant of the current distribution within 
the full counting statistics methodology.— Although the 
measurement of higher order cumulants is experimentally 
challenging, shot noise in atomic contacts and molecular 
junctions has been already measured in the small voltage 
(elastic) regime,—"— and there arc ongoing experimen- 
tal efforts to address the inelastic noise signal as well.— 
These progresses made the investigation of effects due to 
electron-phonon (e-ph) interaction on the current noise 
an appealing task from the theoretical point of view, with 
the ultimate goal of allowing for prediction for the noise 



in molecular junctions as reliable as those now available 
for the non-linear conductance. 

Since the noise is technically represented by a two- 
particle non-equilibrium correlation function, its direct 
evaluation poses a significant challenge compared to the 
mean current. For molecular junctions weakly coupled to 
leads, noise calculations based on one-level models have 
been put forward within the rate equation approach . '^^'^'^ 
In the opposite limit, pioneering works based on the non- 
equilibrium Green's functions formalis m^^'^^ adopted a 
mcan-field-like approximation for the noise, thus neglect- 
ing the contributions due to the vertex corrections. A 
very convenient tool to overcome these limitations is the 
full counting statisticsr^i^ since it allows to address the 
noise and other current cumulants, while taking consis- 
tently into account all contributions due to e-ph coupling 
up to a given order in perturbation theory. Simultane- 
ously with other two groups ;^^'^^ we exploited such an 
approach to analyze the transport properties of a simple 
model for molecular junctions consisting of a single res- 
onant level symmetrically coupled to metallic leads and 
weakly interacting with a single phonon mode.— Despite 
its simplicity, this model can be applied to the exper- 
imentally relevant case of junctions formed by a single 
hydrogen/deuterium molecule between platinum leads^i 
and in this case we predicted a significant inelastic con- 
tribution to the current noise.— 

In this paper we go beyond such a simple model and we 
extend our theory for inelastic noised to more complex 
molecular junctions and to atomic wires. In fact, we con- 
sider the case of a junction formed by multiple electronic 
levels with asymmetric coupling to leads, and derive an- 
alytical formulas for the corrections to current and noise 
due to weak electron-phonon coupling. We express our 
result in terms of universal temperature- and voltage- 
dependent functions and junction-dependent prefactors. 



These expressions, when supplemented with ab-initio cal- 
culations to estimate m.icroscopically the prefactors char- 
acterizing a given junction, can serve as a basis to make 
realistic predictions for the current noise in a relevant 
class of molecular and atomic-size junctions. In this re- 
spect, our work can be viewed as a direct extension of 
the corresponding lowest order expansion scheme devel- 
oped for the non-linear conductance J^^— In addition, we 
identify the contributions to noise due to the vertex cor- 
rections and show that, even in limit of weak e-ph cou- 
pling, they need to be taken into account in order to ob- 
tain accurate results and to comply with the fluctuation- 
dissipation theorem. 

The paper is organized as follows. After a brief de- 
scription of the model of a multi-level junction coupled 
to leads and weakly interacting with a number of vibronic 
modes (phonons) in Sec. |TT1 we introduce the methodol- 
ogy of the noise calculation via extended Keldysh Green's 
functions in Sec. IIIIl In Sec. IIVI we consider the case of 
no interactions and recover the standard results for the 
elastic current and noise. Our original contribution is 
presented in Sees. |V] and IVIl where we discuss the cor- 
rections to the current and noise due to the Hartrce and 
the Fock diagrams, respectively. In subsection IVICI we 
then illustrate our theory by a simple example of indepen- 
dent electronic levels which are coupled only via the e-ph 
interaction. Finally, we conclude and give an outlook of 
open issues and possible extensions of the present work 
in Sec. IVIII More technical parts of the text are deferred 
to 5 appendices. In addition, we make use of the Elec- 
tronic Physics Auxiliary Publication Service (EPAPS) to 
supplement the paper with a Mathematica notebook with 
full expressions for the lowest order corrections to current 
and noise due to e-ph coupling. This file is intended to 
be of use for interested readers in order to access directly 
our results without need of retyping cumbersome formu- 
las from the main text, and it also extends the results of 
subsection IVIBI to the case of finite temperature. 



tonian 



He — H(i + Hph + Hpph 
Hp]i = ^ fujiblbi 



(2a) 
(2b) 

(2c) 



Hepl. = E E ^^'^l'^^ (^1 + ^^)' (2d) 

where di {dj) and b( (6|) are the electron and phonon 
annihilation (creation) operators, respectively; H^ is the 
single-particle effective Hamiltonian of the electrons mov- 
ing in a static arrangement of atomic nuclei, i/ph is the 
Hamiltonian of free uncoupled phonons, -ffgph is the e- 
ph coupling within the harmonic approximation, and 
M^ = {Mff} is the e-ph coupling matrix for the ^-th 
phonon mode. Here, boldface notation stands for matri- 
ces in the system electronic space. The leads and tunnel- 
ing Hamiltonians are given by 



Hl,R = E ^aMci^k^a,k, (3) 

k,a=L,R 

Ht= Y. J2(Kk£lA+h.c.). (4) 



k,a—L,R i 



The states in the leads are occupied according to the 
Fermi distributions /«(£) = /(e — /-'a); with /(e) = (1 -I- 



l/ksT the inverse temperature, and Ha the 



chemical potential of Icad-a. The applied bias voltage is 
eV = fiL - Mfl- 



III. METHODS 

The generalized Keldysh Green's function 
technique 



II. MODEL 



The system we consider can be schematically rep- 
resented as a central device region (representing the 
molecule or the atomic-wire) which is tunnel-coupled to 
non-interacting metallic leads 



H — He + Hl,b. + Hj 



(1) 



Neglecting for simplicity the spin degree of freedomii the 
central region can be described by the following Hamil- 



To calculate the average current and the zero- 
frequency noise in the stationary regime, we employ 
the generalized non-equilibrium Keldysh Green's func- 
tion technique)^ In this approach, one is interested in 
finding the cumulant generating function <S(A), which in 
the case of charge transport is defined as 



-5(A) 



Y.e^'^PtoiN) 



(5) 



N 



where Ptg {N) is the probability for N charges to be trans- 
ferred through the system during the measuring time to 
and A is a continuous parameter known as counting field. 
Given S{X), the cumulants {{6N'')) of the charge transfer 
distribution can be straightforwardly calculated accord- 
ing to the prescription 



{{SN')) = 



Qk 



d{iX) 



^S{\) 



(6) 



A=0 



Under the assumption that the measuring time to is much 
longer than any correlation time of the system (ig — >■ oo), 
the first two cumulants of Pto{N) are directly propor- 
tional to the average current through the system / and 
to the zero-frequency current noise S, 



I = e 



Mm 

to ' 



S = 



to '' 



(7) 



which are the quantities wc arc primarily interested in. 

The key idea for calculating the cumulant generating 
function for transport through a quantum system is to 
modify the Hamiltonian by introducing a time-dependent 
counting field A(t) and to relate iS(A) to the Keldysh 
Green's function of the system in the presence of A(t), i.e. 
to G'l{t,t') = -ih-'^{Tcd^(t)d]{t'))x, where the evolution 
is due to the modified Hamiltonian<^i^ One way to ac- 
complish this is to add a time-dependent phase A(t)/2 to 
the tunneling matrix elements VI j,, 



Hj 



m 






,e 



-«A(t)/2;;t 



c'L,kdj + ykAk'^j + h.c. 



with A(i) = \9{t)d{to — t) on the forward branch of the 
Keldysh contour and X{t) — —Xd{t)9{to —t) on the back- 
ward one,— where 0{x) is the Heaviside step- function. 

Here we extend the result derived by Gogolin and 
Komnik^ for the Anderson model to the case in which 
the central region has several electronic states. General- 
izing the derivation of Ref. |49| to a multilevel system, we 
obtain the following expression for the derivative of the 
cumulant generating function 



dS{\) 
dX 



= t,j ^TrK{K{e)Gx{e)] 



(8) 



where Ga represents the Keldysh-Green's function of the 
system in Keldysh space 



Gx{e) 



G+-(£) G++(e) 



(9) 



Sy = d'Sir/dX, with Sy the self-energy due to the mod- 
ified tunneling Hamiltonian H^ 



St(£) 



iTLe 



-'^'l-.fL{e)]+iTR[l-fR{e)] 



-iTLe^^fLis)-iTRfRie) 



(10) 



r 



and Ttk stands for the trace over the electronic de- 
grees of freedom and the Keldysh space, i.e. Tik{0} = 
Tr{0 -I- 0++}, with Tr{---} being the trace in the 
system electronic space. The check sign " indicates ma- 
trices in the Keldysh space and the superscripts —/+ cor- 
respond to the forward/backward branch of the Keldysh- 
contour. Note that in Eq. (fTO]) we have used the following 
sign convention for the elements of the Kcldysh-matrix 
for the self-energy S 









(11) 



Finally, r^J(e) = 2TrJ2^Vl,^Vl^S{e - ek,a) is the level 
broadening due to the coupling to the lead a. 

According to Eq. ([5]) , the problem of evaluating current 
and noise (as well as any higher order cumulant of the 
charge transfer distribution) is reduced to the calculation 
of the system single-particle Green's function Ga. The 
latter can be obtained from the solution of the Dyson 
equation 

Ga(£) = gx{s) + gA(e)Seph(e)GA(£), (12) 

where Seph is the self-energy solely due to the e-ph cou- 
pling, and gA is the free Green's function of the system 
in the presence of the leads and of the counting field but 
without the e-ph interaction gA = (g^ ^ — '^t)~^ , with 



gd(e) 



el - hd 





d 


el + hd 



(13) 



the Green's function of the isolated dot. It is important 
to notice that Seph, depending on the Green's function 
of the system, is itself a function of the counting field A 
(see Sec. [mC]). 

Finally, we remark in passing that for A 7^ it is 
g^~ -|- g^"^ ^ g^"*" 4- gj~, i.e. in the presence of the 
counting field A, the four Keldysh Green's functions are 
all independent. 



B. Current and Noise 

Although Eq. ([5]) gives access to all cumulants of the 
charge transfer distribution through the system, in this 
work we will focus only on the study of the average cur- 
rent / and the zero frequency noise S, which are the 
quantities most easily accessible from the experimental 
point of view. 

The average current is directly obtained from Eq. (|S]) 
by setting A = 



^/^^"^^^^^("^^^^^^^^-o 



(14) 



while the noise is given by 



S 



— — Trx {S^Ga + S^GaS^Ga}^^q 



2'Kh 



+ ''7 2^TrK{S^GAS^phGA};^^o: 



(15) 



_Q 




FIG. 1. Diagrammatic representations of S// and Sf- Here, 
the plain line stands for the free electronic Green's function 
gA, and the wiggly line for the phononic one df. Finally, the 
dot represents the e-ph coupling constant Mf . 



with A/f = [ef^'^' — 1)~^ the thermal expectation value 
of the mode occupation. The proper inclusion of pos- 
sible heating effects on noise, involving non-equilibrium 
phonon occupation and its potential back-action on the 
electronic transport, is beyond the scope of this publi- 
cation; some of the involved issues are discussed in the 
concluding Sec. IVIIl 

In Eq. P^ we introduced the generalized electronic 
density n^ on the two branches of the Keldysh contour 
[v = ±) in the presence of the counting field 



where, we have used the identity d\d = — 6(9a6 ^)6 
together with the Dyson equation G^^ = g^^^ — 'I^t{^) — 
£eph(A). It turns out that the first term of Eq. ([T5|) cor- 
responds exactly to Eq. (30) of Ref. |50, which gives the 
expression for the noise within a mean-field approxima- 
tion (see Appendix \^. For this reason, we identify 



^(mf) 



' I ^Ttk {S^G. + S^G.S^G,},^^ 



(16) 

as the mean-field contribution to noise. The remaining 
term constitutes the vertex correction 



^(vc) 



As we will discuss in detail in the following, the vertex 
correction S''^'^' can give a significant contribution to the 
total noise, comparable to the mean-field part and thus, 
contrary to what was done in some pioneering works j^i*^ 
it cannot be omitted even in the limit of weak interaction. 
Moreover, neglecting S**-^"^^ generally leads to violation of 
the fluctuation-dissipation theorem, see Appendix [D] 



C. Weak electron-phonon coupling 



In order to make use of Eqs. ([T4| . (|T5|) . we need to de- 
termine the full Green's function G^ . Being interested in 
the experimentally relevant limit of weak electron-phonon 
coupling, we truncate the Dyson equation at the lowest 
(second) order in the e-ph coupling 



Ga ~ gA + gA^gpi^gA 



(18) 



.(2) 



where S^ j^ — S// 4- Sj? is the Hartree-Fock self-energy, 
depicted diagrammatically in Fig. [1] with 






M,Tr{n^Madr(e = «) (19) 



Sr(£) = *E f^dr{e~e')MegT{e')Me, (20) 

and 77,77 = ±. Above, (Cp{e) represent the free thermal- 
ized phonon Green's functions of the ^-th phonon mode 

d±±(e) = E [ - nr{2Nt + l)5{e + stu^,,) ± V ''-—] 

S — ± 

dj-^{e) = -27ri[(7V£ + l)5{e ± nwi} + Ne5{e T h^t)] 



iK= lim -ingr(t-t') 



lim —i 



de 
— ( 

27r 



-ie{t — t'^jhw 



gr(£)- 



(21) 



Note that, on the two branches of the Keldysh contour, 
the electronic density nj^ is given by different limits t ^ t' 
of the corresponding Green's functions g'^^{t — t'). As a 
consequence, even if \{t) = ±A on the forward/backward 
branch, nj 7^ n^^ for any finite value of A. On the other 
hand, at A = one gets 



^A=o ^ "A=0 



i/|ti(e). (22) 



where n^ is the electronic density in the noninteracting 
case. 

Plugging Eqs. P^ . (^0)) into Eq. ([T5|) . we can rewrite 
Eqs. (HI]), (US]) as / = /oi + /eph and 5* = 5oi + 5eph, where 



Ic-\ = ie 



de 
2nh 



TrK{S^gA},^o, 



(23) 



^ci = e^/ ^Trx {S^gA + S^gAS^gAJ^^o , (24) 
are the elastic current and noise, and 

^eph = If + Ih, 5'eph ^ Sp + Sh 

give the respective corrections due to electron-phonon 
coupling, with 



/ 



H{F) = «e 



r de 



and Sh{f) = S^h{F) + ^h{F) ' where 



c(mf) _ p2 



+ 2 S^gAS^gAS^(^)gA| 



(26a) 



A=0 



is the mean-field contribution and 






/5^^^--{^^^^^-(-)^4a.o- (^^^^ 



the vertex correction. The previous equations can be 
schematically expressed by the diagrams of Fig. [5] 



/nl = 



Q(mf) 



= g>., X =t'j. 



^,..— — ««v. 








ji'-— 


"'N^ >) 


b — ^ " 




o 




'S'ol — 


r 


X ) 




. vy. 


A=0 




. v_ 


-^ V^J 


A=0 


/h(F) 


= 


c 


>■ 








. v_ 


.y^. 


A=0 


/^ + 2( 


^ 




c(™) 


/"^ 


V_^ 




%- 


-^ - 


A=0 




. V_^ 


-^ J 



= s;,, n = Sf 



= S' 



FIG. 2. Diagrammatic representations of Eqs. (|23p - (|26p . Tiie 
plain line stands for the electronic Green's function g^; sin- 
gle and doubled crosses stand for S^ and S^, respectively. 
The box represents S^ff^) and, finally, the crossed box stands 
for the derivative of the e-ph self-energy with respect to the 
counting field E'^^^j. 



nance energy £res is far away from the Fermi energy 
(|£ros -Ef\^T, eV, ksT, and Hujq). 

Within the eWBL approximation, the integration over 
energy of functions with compact support can be per- 
formed analytically, resulting in explicit results for the 
mean current and the noise as functions of the applied 
bias voltage and other system parameters. It should be 
noted however, that approximation (ii) potentially leads 
to problems for integrals over infinite range and, in this 
case, it might be necessary to lift it. Specifically, this 
happens in the calculation of the electron density enter- 
ing the Hartree term, see Appendix [B] and in the eval- 
uation of the real parts of the retarded/advanced Fock 
self-energy via Kramers-Kronig relations, see Appendix 

o 



IV. ELASTIC CURRENT AND NOISE 



We note that truncating the Dyson equation to the 
lowest order in e-ph coupling, O(M^), preserves charge 
conservation in that order. This implies that both mean 
current and zero-frequency noise are constant through- 
out the whole circuit, in particular at both junctions be- 
tween the device and the leads. Potential violations of 
charge conservation can only occur in the next order, i.e. 
©(M^), and can thus be safely neglected for any practical 
purposes in the considered weak coupling limit. 

Finally, wc observe that to the lowest order in the e- 
ph coupling, Jeph and Seph are simply given by a lin- 
ear superposition of contributions coming from different 
phonon modes. As a consequence, we can restrict our- 
selves to the case of coupling to a single phonon mode 
with frequency uiq, occupation A/q = [e^^'^° — 1)^^ and 
coupling matrix M. 



D. The extended Avide band limit 

The corrections to current ([25]) and noise (pS)) due 
to the e-ph coupling involve energy integrals which can 
be evaluated in general only numerically. Analytical 
progress can still be made if one assumes the electronic 
structure to be slowly changing over few multiples of 
a typical phonon energy around the Fermi level Ep 
and approximate (i) the level broadening r^ and (ii) 
the non-interacting retarded/advanced Green's function 
gi~ia) -vvith their values at the Fermi energyi^Ti^SiS^ 



r„(£)«r„(SF) = r„ 



^r{a) 



{s)^g^^'^>{E. 



g 



r{a) 



where we took g''^") = gA=o ~ Sa=o ^^ ^^^ defini- 
tion of g'^'-"-'. This approximation, which we call "ex- 
tended wide-band limit" (eWBL), is reasonable for sys- 
tems where either the broadening due to tunneling is 
large (F ^ el/, fc^T, and hojo), or the closest rcso- 



For sake of completeness, before discussing the correc- 
tions to / and S due to the e-ph coupling, we consider 
briefly the results for the elastic current and noise. 

In the eWBL approximation, the elastic current is sim- 
ply proportional to the voltage 



/el = -Tr{T}eF, 



with T = FLg'Tflg", while the noise is given by 



'5'ci — -r 



-TrjT^} + Tr{T(l - T)}U{eV) 



(27) 



(28) 



where we have introduced the function U{x) = 
xcoth{i3x/2). 

The eigenvalues of the matrix T give the "PIN-code" 
of transmission eigen-channels of the molecule connected 
to leads (without e-ph interaction), and Eqs. (|27p . 
(PS)) arc indeed equivalent^ to the standard results for 
current and noise in a non-interacting system derived 
within the scattering theory '^^'^^ However, T is not 
equal to the matrix product tt^ of the transmission 
amplitudes t of the scattering theory (it cannot be 
as T is in general non-hermitian, for example). The 
two matrices are related though by a similarity (non- 
unitary) transformation, which among others ensures 
Tr{T} = Trjttt},^ With this caveat in mind, for sake 
of simplicity we will nevertheless call T the transmission 
matrix in the rest of this paper. The construction of 
the scattering eigenstates within the NGF formalism is 
described in detail in Ref. [SJ. 

We now turn our attention to the corrections to the 
current and noise induced by the e-ph interaction. In 
order to make the discussion as clear as possible, we will 
consider the contributions coming from the Hartree and 
the Fock diagrams separately. 



-,(mf) _ 
^H — 



+ 2 



q(vc) _ 

Su — 



-I A=0 



-I A=0 



FIG. 3. Diagrammatic representations of S]J" and S^'^ . 
Here, plain and wiggly lines stand for the electronic and 
phononic Green's functions, respectively. The single cross 
stands for S^ and the doubled one for Sy. Finally, the dot 
represents the e-ph coupling constant M. 



V. CORRECTIONS DUE TO THE HARTREE 
DIAGRAM 

A. Current 

We start by considering the contributions to the cur- 
rent Ih coming from the Hartree diagram. After inte- 
grating Eq. (j25p in the eWBL approximation we obtain 



Ih = jTr{T(,^^^)}ey, 



with 



-,(qci) 

■H 



2Tr{neM} 

huQ 



TL{g-MAR + h.c.) 



(29) 



(30) 



with He the noninteracting electron density [ Eq. ([22]) ] 
and A„ = g'Tag". The correction Ih is therefore a 
smooth function of the voltage with no features at the 
phonon emission threshold. For this reason Ih has been 
often discarded in previous works on the effects of e-ph 
interaction on the current.— "—'^'^i^'^ 

It should be noticed, however, that Ih is generally non- 
linear in eV, since n^ can be a (smooth) function of the 
applied bias voltage. Such a voltage dependence is nev- 
ertheless rather weak in the eWBL (see Appendix IBJ). 
and in such a case it is possible to interpret Eq. (P^ as a 
quasi-elastic correction to an effective transmission ma- 
trix T = T + T^^ , i.e. because of the e-ph coupling, 
the current is not proportional to the bare transmission 
cocfBcient Tr{T} but rather to Tr{t}. 



B. Noise 

The mean-field contributions and the vertex correction 
to noise due to the Hartree diagram can be schematically 
represented by the diagrams in Fig. [3l which are the re- 
sult of inserting the Hartree self-energy from Fig. [T]into 

appropriate diagrams in Fig. [2 In the usual eWBL, Sj^ 
takes the simple form 






Tr{(l - 2T)T^'='^}C/(el/) + -Tr{TT^°'^}. 

(31) 



Analogously to the current Ih, this contribution has a 
simple interpretation in terms of the renormalization of 
the transmission matrix introduced above T ^ T = T + 
r|n(qo ) ^ Ti^jg (2a,n be seen easily, as Eq. pip corresponds 
exactly to the contribution of order M^ to the elastic 
shot-noise of a system with transmission matrix T: 

^Tr{T2} + Tr{t(l - f)}U{eV) = ^S,i+ 

+ Tr{(l - 2T)T^""}C/(eF) + |Tr{TT<^°"} + ©(M*), 

where ^ei is given in Eq. (^5)) . 

Making use of the cyclic invariance of the trace, the 
vertex correction S]^ can be rewritten as 



S 



(vc) 
H 



e^/h 



2i 
8 



(Tr{MnV})'-(Tr{Mn'_})' 



Re [Tr{Mn'_}] Im [Tr{Mn'_}] 



(32) 



where we have used the fact that [n'_] = — n^, with 
n'^ = {dn'^/dX)\=o. Performing the integrals over energy 
in the usual eWBL approximation^ one obtains 



Re [Tr{Mn'_}] = --Tr{rLAflMg" + h.c.jeV, 
Im[Tr{Mn'_}] = - 



^Trlr^g^MAi - h.c.}^ 



TrlFig^MAfl - h.c.jUieV) 
Ti-{TLiALMAR-ARMAR + h.c:)} 



X y-^-UieV) 



.(mf) 



with gjj = Reg^ Contrary to S]^ ' , Eq. ^ has no 
simple interpretation in terms of an effective transmission 
coefficient and it represents a distinctive contribution to 
noise coming from the Hartree term. From the physical 
point of view, it stems from the coupling of occupations 
of the electronic levels with the current fluctuations 4^i^ 
We note, however, that Im[Tr{Mn'_}] = in the case 
of a system with a single electronic level symmetrically 
coupled to leads. Therefore in this particular case the 
correction to noise induced by the Hartree term is given 

by Sj^ alone. 



VI. CORRECTIONS DUE TO THE FOCK 
DIAGRAM 

A. Current 

We now turn our attention to the corrections to current 
induced by the Fock diagram. Integrating Eq. (PS|) in the 



usual cWBL approximation, wc obtain 



Tr{T^^°'^}eT/ + TriT^^'^jgieV) 



2AAoTr{T^^"'' 



:^"^"}eV^ + Tr{T^ 



e/h 



where 



T^^°"=rL(g'^MgJ,MA^ + /i.c.) 



(33) 



(34a) 



-i(inol) 



= TLg'[MARM--{MAMs-rR-h.c) 



(34b) 

^(asym) ^ ^^^r [^(^^ _ A^jjMg'^rfl + h.c] g'' (34c) 

with A = Al + Afl; = i(g'' — g°) the spectral density. All 
the involved quantities depend only on the properties of 
the system at the Fermi level and can be determined by 
ab-initio calculations^^"— i^^ 

The voltage dependence of Ip is carried by the func- 
tions 

g{eV) = ^[U{eV-nuJo)-U{eV+nuJo) + 2eV], (35) 

and 



HeV) = ]^j de[{np{e) - np{e + eV)) 

X ne'lnpie' -huo) -np{e' + hwo)}ie)].. 



(36) 



where He' {/(£')}(£) = (I/tt) P/de'/(e')/(e' - e) is 
the Hilbert transform. Eq. (j33p is in agreement with 
the result of Viljas et al.^ while a term oc (1 + 

2A/'o)Tr{T^^'''^}el/ is missing in Refs. HI and IM Such 
a discrepancy is further discussed in Appendix [C] 

The functions g{eV) and h{eV) give contributions to 
dl/dV which are even/odd in bias, respectively (see 
Fig. S]). The term proportional to h{eV) vanishes in the 
case of symmetric coupling to the leads, and it is typ- 
ically much smaller than the contribution proportional 
to g{eV), even for asymmetric junctionsi^— Moreover, 
experimentally measured conductance curves are usu- 
ally very weakly asymmetric under reversal of V and 
at present it is unclear if the asymmetry is caused by 
phonons or by other effects. 

At low temperature {ksT <C hujQ, A/q ~ 0), the main 
contribution to Ip is therefore given by the first two 
terms of Eq. ([55)) alone. The first of these terms, lin- 
ear in eV, is a quasi-elastic correction that, similarly 
to Ih, contributes to an effective transmission matrix 



T + T 



(qol) 



. The second one has instead a threshold 



behavior at the phonon emission energy, see Fig. 21 and it 
is responsible for the jump in the conductance observed 
in lETS and PCS experiments. 

The sign of the conductance step at the phonon emis- 
sion threshold (positive or negative) depends on the coef- 
ficient Tr{T)^"'' }, and it has been discussed in detail in 
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FIG. 4. Upper panel: Plots of the dependencies of the func- 
tions g{eV) and h{eV) on the applied bias voltage (Eqs. p5|) 
and (|36p . respectively). Lower panel: Same as above, but for 
the derivatives dg/dV and dh/dV . In both panels fcsT = 
&Jo/30. 



Refs. Il8l - l20i As a rule of thumb, in the case of a molec- 
ular junction with low (high) bare transmission Tr{T}, 
inelastic e-ph scattering results in an increase (decrease) 
of the conductance above the phonon emission threshold. 
In the case of a system with a single electronic level 
symmetrically coupled to the leads via T^ = F^j = F, 



i(inel) 



reduces to T 



(inel) 



(A'/27-2/r2)(i - 2r), where 
T = F^|i?''P is the transmission coefficient. In this case, 
the crossover from an increase to a decrease in the con- 
ductance is predicted to occur at T == 1/2 J^i^ This be- 
havior has been explored and confirmed experimentally 
in Ref . i. 



B. Noise 



Wc finally address the corrections to noise due to the 



Fock diagram Sp = Sp ' - 



■Sp , which are schematically 



represented by the diagrams in Fig. \5\ 

After lengthy but straightforward calculations, inte- 
gration over energy in the usual eWBL approximation 
leads to analytic results for Sp as functions of the 

applied bias voltage. The final expressions are, however, 
rather cumbersome and, for simplicity, we consider here 
only the limit of zero temperature T = 0. The complete 
expressions for Sp at finite temperature are given in 

the supplementary material;^ while the limit eV ^- is 
discussed in Appendix |D] in relation to the fluctuation- 
dissipation theorem. 

In the limit of zero temperature, we obtain 



S 



(mf) 



e2//i 



= Tr{(l-2T)T^^"'^}|eF| 

+ Tr{(l - 2T)T^""''}(|el/| - nwo)e{\eV\ - Ttloo) 



Tr{K^™^^}sign(ey)/i(eT/) 



T=0' 



(37) 
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FIG. 5. Diagrammatic representations of Sp and Sp . 
As in Fig. [S] plain and wiggly lines stand for the electronic 
and phononic Green's functions, respectively. The single cross 
stands for S^ and the doubled one for Sy. Finally, the dot 
represents the e-ph coupling constant M. 



and 



c(vc) 

+ Ti-{K^;''^}sign{eV)h{eV)\^^^, 



(38) 



where T'p"''"'"'' are given in Eqs. (jMal) . (p4b| . 



n(qcl/inGl) 
.(inci) 



+ MAnTLg'-Mg'-Tn + h.c.]. 



(39) 



and 



Kr^ = (l-2T)rL[AflM(Ai-Afl)Mg" + /i.c.] (40a) 



K 



(vc)_ 



M{Ai,TLg'^+g''TLAR)M 
X [ARTL{AL-AR+2ig'R)+h.c.]. 



(40b) 



Finally, 



,s = ±l 



hojQ 



eV 

hujQ 



(41) 



is the zero temperature limit of Eq. ([36|) . The cor- 
rections to noise S 



(nif/vc) 



can then be divided into a 



renormalization of the transmission T 



T + T 



(qol) 



consistently with the sub-threshold correction to the cur- 
rent. Above the phonon emission threshold, \eV\ > huto, 
inelastic processes come into play and their contribution 
to the noise is given both by the second term of Eq. ([37]) 
and by the vertex correction Eq. psp . It is important to 
notice that these two contributions arc in general of the 
same order (sec below the Sec. IVICI for a demonstrative 
example), so that the latter cannot be discarded. 

Experimentally, 4^ is often measured directly by a 
lock-in technique. Such a quantity shows at the phonon 
emission threshold a sharp and distinguishable jump on 
top of a featureless background due to the elastic and 
quasi-elastic contributions. Therefore, we define here the 
inelastic noise signal as the difference of the plateau val- 
ues of the noise derivative just above and below the jump 



dV 



\eV\=hu:o+ckBT 



dS_ 
dV 



\eV\=nwo-ckBT 



(42) 



with c ~ 5 accounting for the finite jump width at finite 
temperatures. At low enough temperatures, terms pro- 
portional to h{eV) give a very small contribution to the 
inelastic noise signal due to the symmetric shape of dyh 
around \eV\ = hujo (for details see Appendix |E]) and we 
can then approximate 



AS" 



-Tr{(l 



2T)Tp'^'> 



Q^"°"}, 



(43) 



symmetric term, which is even in bias, and an anti- 
symmetric one, which contains the Hilbert transform 
h{eV)\r^ and yields an odd contribution. We notice 
that, while h{eV) ^ is a continuous function, its deriva- ^y 
tive shows logarithmic divergencies at eF = zLhajQi^^^ 
These zero-temperature divergencies are, however, an ar- 
tifact of treating the phonons as non-interacting modes, 
and they are regularized either by finite temperature or 
if any broadening of the phonon spectrum is taken into 
account . ^^1^^ This issue, however, goes beyond the scope 
of this work. 

At zero temperature, the symmetric contribution to 
Sp is a piece-wise linear function of eV. At low voltages, 
\eV\ < hojQ, it is given by the first term of Eq. I\'37\i alone. 
Following the same reasoning as for Eq. ([5T|) . this linear 
contribution can be directly interpreted in terms of the 



i.e. at low temperatures AS" carries the structural infor- 
mation about the junction given by the terms with the 
threshold behavior at the phonon emission energy. 



C. Independent electronic levels 

We now consider a toy model for molecular junctions, 
in which we assume the electronic levels to be mutually 
coupled only via the e-ph interaction. In this case, the 
relevant matrices in the system electronic space are given 



L{R)]i] - %rj,L(fl), 



A, + i(r,,L + r,,fl)/2' 



and [M]y = M^, where i,j = 1,...,N, and N is the 
number of electronic levels involved in the transport. Un- 
der the further assumption that each channel is symmet- 
rically coupled to the leads (F/, = Tn = F), the prefac- 

tors TrjK^'"''/'"'^} vanish identically and Eqs. ^, dSH]) 
can be rewritten in a particularly suggestive form in 
terms of the transmission probabilities 71 = r|/(A|+r^) 
of the individual levels 



^(mf) 



N 

2\eV\Y,[^u{l - Ti){l - 27^) + ^7y['n(l - 2%) +7^1 - 27^)] 



j>i 



^ (l--7-)(l-7-) 
Ti % 



} 



N 

+ {\eV\ - huJo)e{\eV\ ~ fu^o) J2 {"^"(^ ~ ^T^)' + 2^7,;,(1 - 271(1 - %) - 27^-(l - 7^))}, 



j>i 



(44a) 



-,(vc) 

Vf 



N 

2{\eV\ - hujo)e{\eV\ - huJo)Y.{2J^^r^{l -7^) +^7^J [71+7^- 27;7;-+ 2y^7;(l -7;)7;-(l -7^)]}, (44b) 



i=l 



j>z 



where we have introduced the dimensionless couphng 
constants 7^ = \AUj\^TiTj/iTiTj). For N = 1, Eqs. g!]) 
reduce directly to the result of Refs. I44l - l46l . The volt- 
age dependence of Sp and Sp is presented in Fig. [S] 
for the case of a systems with only two levels. We no- 
tice that Sp < (see also Ea j44b[) . meaning that the 
vertex corrections correspond to processes that lead to a 
suppression of the noise through the system. Moreover, 
Fig. |6] evidences that the contributions to the noise due to 
the vertex corrections can be of the same order of mag- 
nitude as the mean-field ones, and that they generally 
need to be taken into account in order to make accurate 
predictions for the phonon-assisted current noise. 

In terms of the transmission coefficients of the different 
channels, the inelastic noise signal A5" is given by 

3 ^ 

(45) 



j>i 



with ¥^(71, 7-) - 2[(1 - 7- - Tj? - %{! - 7-) - 7^(1 - 
Tj) — 2-y/7i7j(l — 7i)(l — 75)] . Depending on the values 
of 71 and jij, AS' can be either positive or negative and 
it is in general very sensitive to the parameters of the 
junction as illustrated in Fig. [3 again for the case of a 
system with only two levels. Here we plot AS" as a func- 
tion of the transmission coefficients 7i , T2 for different 
values of the e-ph coupling matrix elements. As gen- 
eral features we notice that AS" is always positive when 
7i , 72 ^ 1 or when they are both close to the full trans- 
mission. Vice versa, AS' is always negative and close 
to maximum suppression for 7i,72 ~ 0.5. Interestingly, 
the characteristics of AS' depend strongly on the rela- 
tive strength of the different matrix elements Afy, and 
therefore the inelastic noise signal might provide a tool 
to extract important information on the e-ph coupling. 



VII. CONCLUSIONS AND OUTLOOK 

In conclusion, in this work we have studied the correc- 
tions due to weak electron-phonon coupling to the aver- 
age current and the zero-frequency noise in a molecular 
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FIG. 6. (Color online) Voltage dependence of Sp (full 
lines) and Sp (dashed lines) for the case of the toy model 



(Eqs. (|44ap and ()44b|) . respectively) for different values of 
the transmission coefficients and of the e-ph coupling ma- 
trix elements. Upper panel: 7i = 0.4, T2 = 0.2, lower 
panel 7i = 0.9, Tz = 0.8. In both panels ksT — 0, and 
Fi = Fa = F. 



junction. To address both quantities in a compact and ef- 
ficient way, we employed the generalized Keldysh Green's 
functions technique. Importantly, for the noise we were 
able to identify distinct terms representing the mean-field 
contribution and the vertex corrections, respectively. We 
included in our calculations both the contributions due to 
the Hartree and to the Fock diagrams and, under the as- 
sumption that the densities of states of the system and of 
the leads depend weakly on energy (eWBL), we derived 
analytic expressions for Ih(f) s-^d Sjjrp\ as functions of 
the applied bias voltage at arbitrary temperature. These 
expressions can serve as a basis for ab-initio calculations 
to make realistic predictions for the current noise in an 
experimentally significant class of molecular junctions. 
Finally, we considered a toy model for molecular junc- 
tions to elucidate the sensitivity of the inelastic phonon 
signal to the parameters characterizing the junction. 

Throughout this paper we have assumed the phonon 
mode to be at equilibrium with an external thermal bath, 
i.e. we have taken the occupation Afo to be fixed accord- 
ing to the Bose-Einstcin distribution Afo = (e^'*"" — 1)~^. 
Such an approximation is strictly consistent with the 
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FIG. 7. (Color online) Contour plots of AS' at zero tem- 
perature as a function of 7i and T2 for different values of 
the e-ph coupling matrix elements Mij. a) Afn ^ and 
M22 = M12 = 0; b) Mil / 0, M22 = O.lMii and M12 = 0; 
c) Mil / 0, M22 = Mil and M12 = ; d) M12 / and 
Mil = M22 = 0. In all panels, Ti = r2 = T. 



lowest-order perturbation theory in the e-ph coupling 
when we implicitly assume strong thermalization of the 
phonon mode. However, it turns out in practice that of- 
ten heating effects cannot be disregarded and that they 
influence in turn the non-linear conductance] ^^i^^i^'^i'^^ 

From the theoretical point of view, the problem of non- 
equilibrium phonon heating can be addressed by extend- 
ing the system Hamiltonian to include the coupling of the 
molecular phonon to other degrees of freedom (typically 
bulk phonons in the leads) . The value of the correspond- 
ing coupling constants can ultimately be obtained from 
ab-initio calculations, which allow to asses the influence 
of the environment from a microscopical descriptioni^ 

In the case of zero counting field A = 0, the non- 
equilibrium phonon occupation for weak e-ph coupling 
can be equivalcntly obtained either by a full non- 
equilibrium calculation evaluating the phonon Green's 
functio n^^'^^1^^1^^ or by solving a master equation de- 
scribing the heating of the devicC j^^'^'^i'^^ which can be 
viewed as a kinetic-equation-like approximation to the 
full non-equilibrium Green's functions studies. Knowing 
the non-equilibrium phonon occupation allows to take 
consistently into account the effects of phonon heating 
in the non-linear conductance. In our pilot study^ we 
used such an ingredient also to phenomenologically in- 
clude heating effects in the noise through a single leveli^ 

However, for a fully microscopical calculation of the 
noise, the situation is considerably more complicated be- 
cause at finite counting field A ^ 0, heating effects cannot 
be expressed solely in terms of the non-equilibrium occu- 
pation of the phonon mode. In fact, to include phonon- 



heating effects in the generalized Keldysh GF technique 
one has to solve the Dyson equation for the phonon 
Keldysh Green's function Dx ~ d + dIlxID\, with the 
polarization operator H^ being given in the lowest order 
by the electron-hole bubbler^i^Si^ see Fig. [8] Note that 
JI\ is explicitly A-dependent via the electronic Green's 
functions and so is consequently also the dressed phonon 
Green's function D\. At A ^ the four Keldysh compo- 
nents of D\ are all independent and therefore, even in the 
kinetic limit (phonon line- width neglected) , it is not pos- 
sible to express the effect of heating just in terms of a sin- 
gle non-equilibrium occupation. It is important to notice 
that substituting D\ for the free phonon Green's func- 
tion in the expressions for the e-ph self-energies S^(p') 
generates extra (additive) contributions to the vertex cor- 
rections 5']^'^ and Sp . These contributions, which are 
related to the influence of phonon fluctuations on the 
electronic transport ( "feedback" ) , are not included in our 
previous phenomenological treatment of heating effects 
on noiso^ and they could possibly account for the dis- 
crepancy between our result and an unpublished one by 
Jouravlev and coworkers?^ which predicts the noise to 
grow with voltage above the phonon emission threshold 
as 5^° ■ -^ ~ T/'*, in contrast to the quadratic behav- 



eph 

ior of Ref. M, S^^^^' ^^^ ~ V^. The idea that phonon 
heating effects could be responsible for a nonlinear volt- 
age dependence of dSeph/dV is further corroborated by 
a recent work by Urban et al.^ which, however, predicts 



S. 



(Rof. 64) 
eph 



V . An independent calculation is therefore 



required to settle this issue. Careful inclusion of phonon 
heating effects into the noise calculations then certainly 
represents a relevant extension of our studies, further- 
more urged by the relevance of heating in several lETS 
and PCS experiments. 

Very recently, a lot of interest has been paid to the 
study of current-induced excitations of local spin degrees 
of freedom in spin-dependent lETS set-ups J^^—i^"— 
Several of these calculation a^^i^^'^^ rely on a perturba- 
tive app roach analogous to the lowest order expansion of 
Ref. [13, also used in this paper. So far, those studies have 
dealt exclusively with the non-linear conductance and the 
study of current noise in those spin systems would be a 
most natural next step. Our method can be straightfor- 
wardly extended in this direction, as long as the occupa- 
tion of spin states is described in a phenomenological way 
via the master equation^ (or just by thermal distribu- 
tion). However, possible further extensions to account for 
fluctuations of a non-equilibrated spin remain, even con- 
ceptually, an open question, because of the anharmonic 
nature of the free spin. Furthermore, the applicability 
of the lowest order expansion itself for the description of 
the spin-dependent lETS experiments seems to be cur- 
rently under debate and renormalized perturbation the- 
ories might be necessary for a proper description of ob- 
served phenomena. Addressing these problems in the 
noise context constitutes an interesting future research 
direction. 
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FIG. 8. Diagram corresponding to the polarization bubble for 
the phonon IIa- The plain lines and the dots represent the 
electronic Green's functions gx and e-ph coupling constant 
M, respectively. 



Finally, the calculation of arbitrary cumulants based 
on the generalized Keldysh GF technique can be imple- 
mented numericallji^iH to address the cases of struc- 
tured tunneling density of states and/or stronger e-ph 
coupling, which go beyond our analytical treatment. On 
the other hand, for the case of weak coupling addressed 
in this work, such numerical methods will face conver- 
gence/efficiency problems due to very sharp phonon line- 
shapes and unnecessary self-consistency loops. In this 
respect, when complemented by ab-initio calculations for 
the transport coefficients Tr{T}, Tr{T]^° },..., our ap- 
proach is designed to be a very efficient alternative to the 
full numerics in the limit of weak coupling and slowly 
varying electronic density of states. It uses the realis- 
tic static calculations of the electronic Green's functions, 
phonon modes, and their coupling as input parameters 
and yields reliable results for the dynamical effects in the 
electronic noise. 
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Appendix A: Mean-field contribution to noise 

Using the invariance of the trace under cyclic permu- 
tations, Eq. (|16p can be recast in the following form 



^ = J dsTv{tTL{fLG> - (1 - /l)G<) 

+ TLG>rLG< + rL{G'' -G^)rL{fLG> (ai) 
- (1 - /l)g<) - /i(i - /i)(G'^riG«ri 

+ G'TLG'Ti)}, 

with G^ = Gj±o, G"" = G^Io - Ga=o and G" = [G'-]t. 
Such an expression corresponds exactlj*^ to Eq. (30) of 



Ref. |50|. We stress that their result was obtained by trun- 
cating the S'-matrix expansion by breaking two-particle 
Green's functions into products of one-particle Green's 
functions, see Refs. ISG orl72t"Sec. 13.8 for further details. 
This procedure holds in a mean-field theory, but it misses 
the contributions given by the vertex correction. For this 
reason we identified Eq. (|16|) with the mean-field contri- 
bution to noise. 

It can be furthermore shown that Eq. (|A1[) is equiv- 
alent to Eq. (10) of Ref. |4l| and to the zero-frequency 
limit of Eq. (9) of Ref. |42|, which therefore represent again 
solely the mean-field contribution to the noise. We stress 
however that, as we will discuss in Appendix |D1 ap- 
proximating the noise with the mean-field contribution 
generally leads to violation of the fluctuation-dissipation 
theorem. 



Appendix B: Electronic density 

The electronic density in the system is given by 

ne = ^^|^g-,U^), (Bl) 



where g^Jp ^ S'^^t g''U=o is the lesser Green's func- 
tion without the electron-phonon coupling. It should be 
noticed that the integrand of Eq. (|Bip does not have a 
finite support and therefore in this case integration over 
energy cannot be carried out in the eWBL approxima- 
tion. Instead, the energy dependence of g''(°) has to be 
taken into account while calculating the integral, and 
only subsequently one is allowed to consider the limits 
r 3> eV, hujQ, UbT or jeres — Ep\ ;2> F, eV, Hujq, ksT cor- 
responding to the eWBL. 

As an example we consider here the case of a system 
with a single electronic level symmetrically coupled to 
unstructured leads with constant F^ = F^; = F. In this 
case 



rie. = 



. f de Tihje) -f /fl(£)] 



Assuming zero temperature and symmetric voltage drop 
at the barriers /i^ = —fiR = eV/2 one gets 



1 1 

2^2^ 

1 1 

2 " ~ 



arctan 



eV - 2e, 



2F 



arctan 



eV + 2eo' 



Teo 



TT V F y 47r(F2 + £2)2 



2F 

{eVr + 0{{eVf), 



which shows that in the eWBL He depends very weakly 
on the applied bias voltage. 



Appendix C: Explicit form of the Fock self-energy at 
zero counting field 



At A = 0, the Keldysh components of Sf satisfy the 
identity (S^^ + S++)a=o = (S^+ + S^")a=o. In this 
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case, it is meaningful to introduce the retarded and ad- 
vanced self-energies by SJ, = (S^~ — S]^^)a=o and 

Tip = [^f] • For definiteness, we also introduce the 
notation T^ = SJzq for the lesser and greater compo- 
nents at zero counting field. The latter can be easily 
calculated from Eq. ([20]) giving 



is the origin of the discrepancy between our result for 
If, Eq. ([55)) . and the expression derived by Paulsson 
and coworkers . ^^1^'^ This discrepancy stems from the sub- 
tleties in the use of the eWBL mentioned above, and was 
already pointed out by Viljas et al.r^ whose result agrees 
with ours. 



a=LM 

+ (AAo + 1)A„(£ + ncjo)./a(e + nLOQ)\m 
SF(e)=-* Y. M[A/-oA„(e + ritJo)(l-/a(e + ?i^o)) 

a=L.,R 

+ (AAo + 1)A„(£ - nw„){l - /„(£ - nwa))]M. 

The retarded self-energy can in turn be written in terms 
of the lesser and greater components using the identity 
S^ — T'p = T,p — Tp and Kramers-Kronig relation 
ReSJ,(e) = He' {imS J, (£')}(£)■ This leads to 

ImS''(e) = -iM{(AAo + l)A(e-fi£^o)+A/'oA(e+fiwo)}M 

a=L,R 
- Aa{e - hiJo)fa{£ - /iwo)}M, 

ReS''(e) =M[{^fo + l)gjj(e - fiwo) + A/'ogfl(e + fiwo)]M 
-^ E M[n,'{A^{e')Ue')}{e + nuo) 

a=L.R 

-n,,{A„{e')Ue')}{e-nuo)]M, 

where we have used the identity Hs'{A{e')}{e) = 
-2ne'{g'i{e')}{e) = -2g^(e), with g^^^^ the real (imag- 
inary) part of g''. We point out that the energy 
dependence of g''^") cannot be disregarded while us- 
ing Kramers-Kronig relations, as the Hilbert transform 
He' generally involves integrals over infinitely extended 
range. However, in the limits F 3> eV,huja,kBT or 
£,08— ^_f| 3> r, eV, hajQ, UbT corresponding to the eWBL 
approximation the previous expressions take a simpler 
form 

ImEJ,(£)=-i Y. MA„M[(2AAo + l) 

a=L,R (CI) 

+ .fa{£ + tlUJo) - fa{e - fiiJo)] , 

ReS5,(e) = (2A/-o + l)Mg^M-i ^ MA„M 

a=L.R (C2) 

X ne'iUie' + hiuo) - Us' - hujo)}{e). 

Note that the Hilbert transform now involves only a func- 
tion with finite support. Inserting these expressions into 
Eq. (PSj) , it is easy to show that the first term of Eq. (|C2[) 



Appendix D: Fluctuation-dissipation theorem 

The fluctuation-dissipation theorem relates the noise 
at zero voltage to the linear conductance of the system 
G 



S{V = 0) = -G. 

In the case of the contributions due to the Hartree 
term, it follows form Eqs. ([32]), (jSl]) that S'jj''-' -^ 

at zero voltage, while S^ fulfills the fluctuation- 
dissipation theorem 

SHiV = 0) = s'^'\v = 0)^lGH, 



P 



^(qoi)i 



with Gh = eV/iTr{T^°'^}. 

The situation is, however, different for the Fock term. 
In fact, in the limit eV — >■ 0, both the mcan-fleld contri- 

.(mf) 

'f 
from zero and reduce to 



bution S'^p"' and the vertex corrections Sp are different 



i(mf) 



(y = o) 



e'^/h 



Tr{T^^"'^} 



(qei). 2 u{hLoo) u{hLJo)^-{huoy 



(3 hojo 



fiwn 



X TrjT*^""'^ + r^g'^MAiMg'^}, 



s'^^\v^o)_ ^^,^ „.,,, ^,„a,u{f^or~{nuoy 



e^/h 



-TrlFLg'-MALMg"}- 



hhjQ 



On the other hand, the correction to the linear conduc- 
tance due to Sf is given by 



Gf 



9/;'^ IF J t ^-Lri-Lp j . 



Comparing the previous expressions one can see that 

sp^\v = 0) + s^;''\v = 0) = ^Gf. 

but the mean field contribution alone does not satisfy the 
fiuctuation dissipation theorem Sp (V = 0) ^ 2Gf/ P- 
This clearly shows that in general, even in the limit of 
weak e-ph coupling, vertex corrections must be included 
into the noise calculation in order to obtain consistent 
results. 
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where Kj^™ '^'^'^ are those of Eqs. (|40|) and 



-(mf) 



K^' = riAflFiiAHM (Ai - A;,J g'' + /i.e.], (E2a) 

K^"^> = iM(A;,rLg'^+ ;j.c.)M(ArLg^- /i.e.) - k1"^V2. 

(E2b) 



The Hne shape of S^yj^ is defined by the functions 



h-^{eV) = /i(eV)coth(/3ey/2), 



(E3) 



FIG. 9. Upper panel: Plots of the dependencies of the 
functions /ii(eF) and h-iieV) on the applied bias voltage 
(Eqs. (|E3P and (|E4p . respectively). Lower panel: Same as 
above, but for the derivatives dhi/dV and dh2/dV. In both 
panels ksT = fiwo/30. 



Appendix E: Anti-symmetric contribution to Sf 



with h{eV) given in Eq. ()36p and 



/i2(eV") ^ ds npie + eV)[l - npie + eV)] 



(E4) 



In this appendix we give the complete expression for 
the terms of Sp which are anti-symmetric with re- 

spect to the bias voltage. At finite temperature they are 
given by 

^asym _ TrlK^"')}/!! (eT/) + TrjK*'"' ^ }/i2 (eF), (Ela) 



C'(vc) 
^asym 

eyh 



-(vc) 



Tv{K';''>}hi{eV) + Tr{K)^'>}h2ieV), (Elb) 



(vc)i 



We notice that h2ieV) = {eP)-^dvh{eV), i.e. h2ieV) 
is directly proportional to the derivative of h{eV) but 
it exhibits no divergencies and actually tends to zero in 
the limit T ^- due to the suppression factor 1//3 with 
respect to h{eV). Furthermore, we observe once again 



that SI 



(mf /vc) 
asyin 



give a negligible contribution to the inelas- 



tic noise signal AS" at low temperatures, since dyhi is 
almost symmetric around the phonon emission thresh- 
old, i.e. dyhil, ,,|^j, 
suppressed by low temperature, see Fig. [HI 



^v/il|leV|<ftc.c' ^'^'^ ^^^"2 is 



federica.haupt(a)uni-k onstanz.de| 
tno@karlov.mff.cuni.czl 



■'■ wolfgang.belzig@uni-konstanz.de 

^ G. Cuniberti, G. Fagas, and K. Richter, eds.. Introducing 
Molecular Electronics (Springer, Berlin, 2005). 



^ N. Agrait, A. L. Yeya ti, and J. M. van Ruitenbeek, Phys 
R|ep., 377, 81 (2003)|, ISSN 0370-1573 



^ V. G. Naidyuk and I. K. Yanson, Point-Contact Spec- 
t roscopy (Springer, Berlin, 20015). 



^° T. Frederiksen, K. J. Fra nke, A. Arnau, G. Schulze, J. I 
Pascual, and N. Lorente, |Phys. Rev. B, 78, 233401 (2008) 



,lh 



M. Rahimi and M. Hegg, phys. Rev. B, 79, 081404 (2009) 

^^ C. R. Arroyo, T. Frederiksen, G. Rubio-BoUinger, 

M. Velez, A. Arnau, D. Sanchez -Portal, and N. Agrait. 



phys. Rev. B, 81, 075405 (2010) . 
^^ A. J. Heinrich, J. A. Gup ta, C. P. Lutz, and D. M. Eigler, 



science, 306, 466 (2004) 



R . G. Jaklevic and J. Lambe, Phys. Rev. Lett., 17, 1139 



(:966 



C. F. Hiriibehedin , C. P. Lutz, and A. J. Heinrich, 



TS- 



312, 1021 (2006) 



Science 



B C. S tipe. M. A. Rezaei, and W. Ho, Science, 280, 1732 



Y.-S. Fu, T. Zh ang, S.-H. Ji, X. Chen, X.-C. Ma, J.-F Jia, 



and Q.-K. Xue, phys. Rev. Lett.. 103. 257202 (2009) 
^"^ N. Lor ente and M. Persson, phys. Rev. Lett., 85, 2997 



Agrait. C. Untiedt. G. Rubio-BoUinger, and S. Vieira, 



(2000) 



Phys. Rev. Lett., 88, 216803 (2002) 

■^ R. H. M. Smit, Y. Noat, C. Untiedt, N D. Lang, M. G 

van Hemert, and J. M. van Ruitenbeek, [Nature, 419, 906 



(:i002) 



M. Brandbyge, N. Lorent e, and A.-P, 



Jauho, phys. Rev. Lett., 93, 256601 (2004) 
^^ M. Paulsson, T. Frederiks en, and M. Brandbyge, Phys 



D Djukic, K. S. Thygesen, C. Untiedt, R. H M. Smit 



K . W. Jacobsen, an d J. M. van Ruitenbeek, Phys. Rev 



Rev. B, 72, 201101 (2005) 
" J. K. Viljas, J. G. GuwasT^ F. Pauly, and M. Hafner, phys 



Rev. B, 72, 245415 (2005) 



E |, 71, 161402 (2005)| . 
^ O. Tal, M. Krieger, B. Leerink, and J . M. van Ruitenbeek, 



L. de l a Vega. A. Martir>Rodero. N. A grait. and A. L. 



Phys. Rev. Lett., 100, 196804 (2008) 



Yeyati, phys. Rev. B, 73, 075428 (2006) 



14 



G. C. Solomon, A. Gagliardi, A. Pecchia, T. Fr a,uenheim. 
A. D. Carlo, J. R. Reimers. and N. S. Hush. []L_CheiiL 
P |hys., 124, 094704 (2006| 



spin degeneracy is explicitly tak en into account. 

^^ L. S. L evitov and M. Reznikov, [Phys. Rev. B, 70, 115305 



N. 



7 3, 233418 (2007) 



(2004) 



Sergueev, A. A. Demkov, and H. Guo, IPhys. Rev. B, "*" A. OT Gogolin and A. Komnik, Phys. Rev. B, 73 



(2006) 



T. Frederiksen, M. Paulsson, M. Brandbyge, and A.-P. 



195301 



J mho, Phys. Rev. B, 75, 205413 (2007) 



F. M. Souza. A. P. Jauho, and J. C. Egues, Phys. Rev. B, 



78, 155303 (2008) 



G. Teobaldi, M. Penalba, A. Arnau, N. Lorente, and W. A. 



Hofer, 



Phys. Rev. B, 76, 235407 (2007) 



Y. Me ir and N. S. Wingreen, Phys. Rev. Lett., 68, 2512 



(1992) 



T. Fr ederiksen, N. Lorente, M. Paulsso n, and M. Brand 



b yge, [Phvs. Rev. B. 75. 235441 (2007) , 

M. Paulsson, T. Frederiksen. H. Ucba, N. Lorente, and 



'"^ M. Biittiker 

'^^ M. Pau lsson and M. Brandbyge, Phys. Rev. B, 76, 115117 



Phys. Rev. B. 46, 12485 (1992 



M. Brandbyge, phys Rev. Lett 100 226604 (2008)1 
L S. Kristensen' M. Paulsson. K. S. Thygesen, and K 



TJT- 



(2007) 



J. Fran sson and M. Galperin, Phys. Rev. B, 81, 075311 



Jacobsen, Phys. Rev. B, 79, 235411 (2009; 



W. 



(2010) 



M. Alducin. D. Sanchez-Portal. A. Arn au. and N. Lorente, 
Phys. Rev. Lett., 104, 136101 (2010). 



33 

34 



J . Fransson. H. C. Ma noharan. and A. V. Balatsky, 



3tt., 10, 1600 (2010) 



Nano 



56 

57 



Unlike lie, the evaluation of n^ involves only integrands 
with compact support . 
S. Hershfield, 



Phys. Rev. B, 46, 7061 (1992)| . 



S . Monturct. M. A lducin. and N. Lorente, Phys. Rev. B 



8 I, 085447 (20 10)|. 

K. R. Patton, arXiv:100V.1238yl (2U1U), (unpublished) 

d 



See associated Mathematica notebook for the complete ex- 
pression for Sp at finite temperature. 



O. Entin-Wohlman, Y . Imry, and A. Aharony, Phys. Rev 



B, 80, 035417 (2009) 



Y 



Beenak ker and C. Schonenberger, Physics Today, 56 



3 " (2UU3) 



R. Egg er and A. O. Gogolin, Phys. Rev. B, 77, 113405 



(2008) 



M. Blanter and M. Biittiker, Phys. Rep., 336, 1 (2000) . 

.. V. Nazarov, Ann. Phys., k, «!-!!):] [im); V. V^. ^ 

Nazarov, ed.. Quantum Noise in Mesoscopic Physics ^^ 



M. Engelund. M. Bra, ndbvge. and A. P. Jauho, Phys. Rev 



B. 80, 045427 (2009) 



Y. Asai, Phys. Rev. B, 78, 045434(2008} 



( Springer, Berhn, 2003). ^ 

h. IS. van don Jjrom an d J. M. van Kuitenbeeit. Phvs. Kev 



D. A. R yndvk and G. Guniberti, Phys. Rev. B, 76, 155430 



;tt., 82, 1526 (1999) 



(2007) 



p. Diukic and J. M. van Ruitenbeek, Nano Lett., 6, 789 



(2006) 



"■^ Its generalization to multilevel case is st raig htforward and 
follows exactly the lines of Refs. [13 and [23. 



64 



D. F. Urban, R. Avri ller, and A. Levy Yeyati, Phys. Rev 



B, 82, 121414 (2010) 



Kiguchi. U. lal. S. Wohlthat. 1. Paulv. M. K riegerT 



D. Djukic, J. G. Cuevas, and J . M. van Ruitenbeek, Phys 

r[ ■ 



ev. Lett., 101, 046801 (2008' 



"° O. N. Jouravlev, Noise and Spin in Nano structures, Ph.D 
thesis, TU Delft, Nethe rlands (2005). 



66 



M. van Ruitenbeek, private communication. 



A. Mitra. 1. A leiner, and A. J. Miilis. [Phvs. Key. B. 69 
2^5302 (2004)[ 



J. Fer nandez-Rossi er. Phvs. Rev. Lett.. 102. 256802 



(2009)1; M. Persson, Vbid., 103. 050801 (2009)| ; J. Fransson, 



■. Wm-. J. KoTK 



Koch and F. von Oppon, Phys. Rev. Lett., 94, 206804 



M. E. Ralkh. and V 



(Jppen, \ibid. 



[Nano Lett., 9, 2414 (2009) 
^^ N. Lorente an d J. -P. Gauyacq, Phys. Rev. Lett., 103, 



9 P7056801 (2005) 
*^ J .-X. Z hu and A. V. Balatsky, [Phys. Rev. B, 67, 165326 



176601 (2009): J.-P. Gauyacq. F. D. Novaes. and 



( j003)[ 
*^ M. Galperin. A. N itzan. and M. A. Ratner, [Phys. Rev. B 



7 1, 075326 (2006) 



N. Lorente, phys. Rev. B, 81, 165423 (2010)[ F. D. Novaes 
N. Lorente, and J.-P. Gauyacq, \.bid., 82, 155401 (2010) . 
®* F. Delgado, J. J. Palacios, and J. Fernandez-Rossier, 
104, 026601 (2010)[ F. Delgado and 



Phys. Rev. Lett 



J. Fernandez-Rossier, arXiv:1006.5608yl (2010) 



S. Levitov and G. B. Lesovik, JETP Lett., 58, 230 



J. Fransson. O. Eriks son, and A. V. Balatsky, Phys. Rev 



(1993); L. S. Levitov, M. W. Lee, and G. B. Lesovik, 
J. Math. Phys., 37, 4845 (1996) 



T. L. Schmidt and A. Komnik, 



009) 



Phys. Rev. B, 80, 041307 



' R~^r iller and A. Levy Yeyati, Phys. Rev. B, 80, 041309 



009) 



■'^ F pHaupt. T. Novot nv. and W. Belzig, [Phys. Rev. Lett. 
l t)3, 136601 (2009)[ ' 



B, 81, 115454 (2010)[ . 
F. Elste and G. Timm, [Phys. Rev. B, 81, 024421 (2010)| . 
Apart from an extra factor of 2 in Eqs. (30), (31) of 
Ref. [5y stemming from their different definition of the zero- 
frequency noise. 

H. Haug and A.-P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors, 2nd ed. (Springer, Berlin, 
2008). 



*^ Based on a spin-less model, our results need to be multi- 
plied by a factor of 2 when compared with works where 



